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Abstract 

In order to check the validity of auxiliary field method in the Nambu-Jona-Lasinio 
model, the one-loop (=quantum) effects of auxiliary fields to the gap equation are consid- 
ered with iV-component fermion models in 4 and 3 dimensions. N is not assumed so large 
but regarded as a loop expansion parameter. To overcome infrared divergences caused by 
the Nambu-Goldstone bosons, an intrinsic fermion mass is assumed. It is shown that the 
loop expansion can be justified by this intrinsic mass whose lower limit is also given. It is 
found that due to quantum effects, chiral symmetry breaking (xSB) is restored in D = 4 
and D = 3 when the four-Fermi coupling is large. However, xSB is enhanced in a small 
coupling region in D = 3. 
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1 Introduction 



The Nambu-Jona-Lasinio (NJL) model^X is, needless to say, one of the most famous 
field theoretical model exhibiting dynamical chiral symmetry breaking (xSB) phenomena. 
Owing to its simplicity, there have been many studies using the NJL models: some of 
recent trends are those which deal with external disturbances such as background gauge 
fields [21 IU 03 or curved spacetimejHj in order to explore the detailed phase structure. 
Originally, calculations for dynamical xSB had been done in a self-consistent manner to 
yield the gap equation but later it was revealed that results are much easily obtained 
in path integral with the use of the auxiliary]?] or the Hubbard-Stratonovich 8 field. 
The recipe which we shall call the auxiliary field method jSj becomes exact when N, the 
number of degrees of freedom of the original dynamical fields, goes to infinity. However, 
the analysis in lower dimensional bosonic models shows that even if N = 1, we can improve 
results toward the true value by taking higher orders in the loop expansion lO . Therefore 
it is desirable to incorporate quantum(=loop) effects of auxiliary fields to the gap equation 
when N remains finite. 

However, there is an obstacle to perform the loop expansion in terms of auxiliary fields 
in the NJL model: due to massless Nambu-Goldstone bosons, members of auxiliary fields, 
infrared divergences are inevitable in higher loop calculations. Kleinert and Bossche pay 
an attention to this infrared regime to conclude that there is no pion in the NJL model[TT]. 
that is, there is no room for the auxiliary fields. Their main ingredient is, however, a chiral 
nonlinear model 12 , an effective theory, so that there need more rigorous and careful 
investigations. Indeed some oppositions to this conclusion have been raised [13] • 

We follow the standard prescription for the effective action formalism by introducing 
sources coupled to bilinear terms of fermions. In order to control the infrared singularity, 
we assume an intrinsic mass of fermions, so called a current quark mass. We do not care 
about the renormalizability so that an ultraviolet cutoff is introduced to define a model. 
We will make no approximation other than the loop expansion. 

There have been attempts to consider 0(1/N) terms in four-Fermion models: some 
calculate the effective potential for the model with a discrete chiral symmetry in order 
to clarify the renormalizability in 4 > D > 2^1] , others check the Nambu-Goldstone 
theorem, Geil-Mann-Oakes-Renner or Goldberger-Treiman relations in the 4-dimensional 
NJL model 1 to show that 0(1/N) effects weaken quark condensation (qq) obtained by the 
tree order ^j. However, there seems no attempt to study the higher loop contribution, 
paying an attention to the infrared regime, to the vacuum condition and to the gap 
equation in terms of auxiliary field method. In §2, we present a general formalism in 
path integral to obtain an effective potential in the NJL model. We work with the N- 
component fermion model in D = 4 and D = 3 16 but as stated above N is merely a 
loop expansion parameter to be kept finite. The next section §3 deals with the vacuum 
condition and then the gap equation up to the one-loop order of the auxiliary fields. It 
is concluded that quantum(=loop) effects restore ^SB in D = 4 while the situation in 
D = 3 is slightly different; x^B is restored in a strong coupling regime but enhanced in a 
weak coupling regime. We also find the lower limit of a current quark mass to ensure the 

1 It is, however, almost trivial to check these relations under 0(1/ N) expansion; since it is well known that 
those are the consequences of the Ward-Takahashi relations(WTs) which are persistent with the loop(=C(l/iV)) 
expansion. 
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loop expansion. The final section is devoted to discussions. 



2 Model and Basic Formalism 



In this section, we develop a general formalism in order to clarify our goal: the NJL Model 
with an intrinsic mass e in 3 as well as 4 dimensions is given as 



Ant 



A 

2N 



1>(x) (P + e) + Ant 

(V>0)V>(x)) 2 + {i>{x)i^{x))' 
(ip(x)tp(x)) 2 + (ip{x)i~f^(x)Y + (VK^H^VK^))' 



(1) 



D = 4 
D = 3 



A.(^ x )r^(x)) 2 ; 



with 



r a = (i,ir 5 ) ■ r 5 



75 , D = 4 

(74,75) , £> = 3 



(2) 



(3) 



where ft = 7^^ and iV-component fermion fields have been introduced and 7^'s are 4x4 
matrices, even in 3 dimension 2 , satisfying 



{l^lv) =26^; n,v = 1,2,3,4,5 . 



Explicitly, 

7m = 



^ 
-a,. 



; /•* = 1 ~ 3 , 74 



1 

1 



75 = 71727374 



a 
-a 



(4) 



The intrinsic mass e, so called a current quark mass, has been assumed to prevent an 
infrared divergence in pion loop integrals. (See the followings.) 
The quantity we should consider is 



Z[J] = Tr T exp 



dtH{t) 



H(t) = I d D - x x 



N 



lP{x) ( 7 • V + £ - J(X) • r) lP{x) - Ant " ^J\x) 



J(x) = J a (x) 



(J(x),J 5 (z)) , D = A 

(J(x),J A (x),J 5 (x)) , D = 3 



(5) 
(6) 

(7) 



In 3 dimension, we need an additional (iV-component) fermion to form a four-component spinor [Trjl H] 

i> = ( j . i> = V> f 73 = ( iPi -^2 ) = ( V4 CT 3 -V'ks ) , 



to be able to realize a chiral symmetry, (when e = 0), 

i/>(x) — ► e ia ^ij}(x) , ip(x) 



ci/375 



ip(x) 



which is, therefore, a global U(2) symmetry, finally broken down to U(l) x U(l) by a mass term. 
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where J(x)'s are c- number sources and T designates the (imaginary) time-ordered prod- 
uct. (We have introduced a J 2 term just for a later notational simplicity 14 .) From this 
we can extract an energy of the ground state by putting J i— ► as well as T i— ► oo. The 
path integral representation reads 



Z[J] = J d[ip]d$]exp I J d D x 
Introducing auxiliary fields in terms of the Gaussian integration, 

I / d[27]exp 



(8) 



S(x) = S a (x) = 



(a(x),vr(x)) 



D = 4 



v (a(x),TT 1 (x),7r 2 (x)) , D = 3 
to eliminate the four- Fermi interaction, we find 

N 



(9) 
(10) 



Z[J] 



d^]d[tp]d[U] exp / d 



2A 



N 



-^{^ + e + {E-J)-r}^ + ^J 2 



(11) 



The fermion integrations yield 



d[U] exp 



-JV 



J d D x^(£ 2 - J 2 ) - Trln{@ + e + (E - J) ■ T} 



j d[E]e W (-NI[U,J}) , 



with 



d D x[ —U Z + ^E-J 



1 



1 



Trln {fi + e + E- T) 



2X A 

where Tr designates spinorial as well as functional trace. Write 



Z[J] 



-NW\J] 



and introduce "classical" fields, 



1 J 5W eq.fO 1 „ eq. 

a^TT = \ <E> = 



(12) 



(13) 



(14) 



(15) 



with expectation values being taken under the expression (|12(l or under the original one 
(|5|). then perform a Legendre transformation with respect to W[J] to obtain the effective 
action, 



r[<f>] = w[j]-~[j 



(16) 



where a shorthand notation, 



A • b) = j d D x A(x)B(x) , 



(17) 



has been employed. When «7's are set to be constants, an effective action becomes the 
effective potential: 



r[<f>] J ^£ st - VTV(cf>) 



(18) 



with V being (D — l)-dimensional volume. 

We calculate W[«7] with the help of the saddle point method 3 : first, find the classical 
solution, Eq, 







SI 



SE(x) 



So 



(Z + J)(x)-trrS(x,x : S ) , 



where S(x,y : Uq) is a fermion propagator under the background fields, 

(P + e + E Q (x)-r)S(x,y:E ) = 8(x-y) . 
Second, expand I around Uq 



1 



3! 



r(3) 



where 



r(n) 



5 n I 



5U n 



So 



(4 n) • 



j d D Xl ■ ■ ■ d D x r 



5 n I 



5£ a ^ {xi)---5E a ™(x n ) 



(19) 

(20) 
(21) 
(22) 

U^( Xl )---U a -(x n ) . (23) 



Third, put (X 1 — Uq) S/y/N, then perform the Gaussian integration with respect to 
U to obtain 



^]='o+^Trln4 2) +0(^) 



where 



2AV 



Z -Z Q ) +- [So- J) -Trln 



+ ^0-^) ; 



(24) 



(25) 



(4 



(2) 



5 2 I 



ab 5U a (x)SE b (y) 



So 



-5{x - y)5 ab + tvr a S(x,y : U )r b S(y,x : S Q ) .(26) 



(2) 

Here the trace is taken only for spinor space. Jq is a matrix in the £ space (2 x 2 in 

(2) 

D = 4, 3 x 3 in D = 3). Jo is the "tree" part while Trln/p is the "one-loop" part of the 
auxiliary fields. Using eq. (|15j) we find 



0(x) = S (x) + 



A 



2N SJ(x) 



Trlnl^) = S (x) + 



gi(g) 
N 



(27) 



3 The result becomes exact when N goes to infinity, which, however, is not the case in this analysis: N is a 
mere expansion parameter that finally be put unity. 
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Note that difference between cf) and Uq is 0(1/N). Inserting ea.l|27jl into the effective 
action (fTTTf) with the use of eqs. (|24|) . (|2"5|) . and (|2"?t|) . we obtain 



r[4>] = 2A^^J -Trln^ + e + ^-T) 

+ ^Tr In (^(z - y) + tirS(x,y : cf>)rS(y,x : 0)J + O(^) • (28) 

By putting J's into constants, the effective potential (fT%|) reads 

V(0) = ^0 2 -^T¥ln(^ + e + 0-r) 

+ 2ivW Tr ln " y) + trr,s(x ' y : 0)r5(y ' x : 0) ) + ■ (29) 

In ea. (|29|) the first two terms are the tree part and the third is the one-loop part of the 
auxiliary fields, whose functional trace should be also taken for the £ space. The vacuum 
is chosen by 



dV_ 



= . (30) 

j=o 



Armed with these, we now proceed to a detailed calculation. 

3 Vacuum and the Gap Equation 

Write the "classical" fields at J = as 

<t>=(m,Z„); «-={ ( ^), D D Zt\ ■ (3D 

and 

4> = (m + e, £ w ) , (32) 
to study the vacuum condition (|3fl|) on the tree part of the effective potential (f29|) , 



iK + ^-f/^ln^W), (33) 



giving 



dVn m M , . f d p 1 

si = ° = T- 4 <'"W(^V7^ < 34 > 
^ 0= * <W,gp i . (35) 



dS v A "7 (2vr)^ p 2 + ^2 • 



6 



(Recall trl = 4 for both D = 4 and D = 3.) Therefore the solution is 

2* = , 



Tfl f 

— = 4(m + e) 



d D p 



1 



(2vr)^ p 2 + (m + e) 2 ' 
whose second relation, called the gap equation in the tree order, reads 



where 



1 

Ad 



/a; — e 



x = 



(m + e) 

7P~ 



<x<1) , 



(36) 
(37) 

(38) 
(39) 



and 



<?i i x ) = 1 — x In 
2 



1 + x 



1 — \/x tan 1 —= 

x 



X D 



XA 



D-2 



(40) 
(41) 



4 D/2-l 7r D/2 r ( jD/2 ) ' 

with T(D/2) being the gamma function; T(2) = l,T(3/2) = y/n/2. In view of en.(l38l 



(42) 



X^O^-e) , 

which implies a trivial fact that mass is e even in the free theory A = 0. 
The one-loop part Vi of the effective potential (|2*U)l reads 



(43) 



(44) 



where tr should be taken for the spinorial space as well as the S. In ea.(|44|) the argument 
of the logarithm is nothing but a two-point function of the auxiliary fields, therefore, we 
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call II a vacuum polarization matrix, 
d D l ( 1 



n , , f d u l 



1 



d D i 



1 



-r, 



x< 



it 



{1 + V/2Y + & (l-p/2) 2 + fo 



( ( -(l 2 - P -) + { m + ef-^ 



V 



2(m + e)ir 



2 P 



2(m + e)vr 

2 



/ -(/ 2 -^) + (- + ^) 2 -^ 
2(m + e)7Ti 
2(m + e)7r 2 



(m + e) 2 + 7T 2 



2{m + e)7Ti 



L> = 4 



-(/ 2 -^)-(m + e) 2 + vr 2 -. 2 



(45) 



27T17T2 



2(m + e)7r 2 

27T17T2 



\ 



By noting 



tr 



2 P 



(m + e) 2 - irf + tt| 



; D = 3 . 



an 



and the tree relation (|35|) . the vacuum condition for XV up to the one- loop is 



(46) 



<9V 



J=0 



(47) 



with V = Vo + Vi , allowing us to choose U n = as a vacuum. The second fact is that in 
the one-loop part of the gap equation, 



dVi 



dm 



(48) 



J=0,S n =0 



we can utilize the tree result (|H7[). 

Therefore in the gap equation we put S n = in the expression (|45|) to find the diagonal 
matrix 



n 



^=0 



ni o 

o n 2 

ni o o 

o n 2 o 
o o n 2 



D = 4 



D = 3 



(49) 
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with 



Hi 

n 2 



d D i 



(2n) D {(I + p/2) 2 + (m + e) 2 } {(/ - p/2) 2 + (m + e) 2 } 



Write 



l 2 - P - 



p 2 _ (I + p/2) 2 + (I -p/2) 
4 " 2 



(50) 



(m + eY 



2 2 
p 



then 



III 

n 2 



+ 



(2vr)^ (/ + p/2) 2 + (m + e) 2 7 (2vr)^ (/ _ p / 2 ) 2 + (m + e) 2 



d D l 





( (m + e) 2 +p 2 /4 \ 


d D l 


\ p 2 /4 } 



{2tt) d | (/ + p/2) 2 + {m + £)2 | |(/ _p/2) 2 + ( m + e ) 2 } ' 
Shifting the momentum (although we are in a cutoff world), we obtain 

: -4 



(51) 



III 

n 2 



i 



+8 



(27r) D P + (m + e) 2 

1/2 



1/2 





/ (m + e) 2 +p 2 /4 \ 




I P 2 /4 / 



{2ir) D [p + p2(i/4 _ t 2) + ( m + £ )2] 



21 2 ' 



(52) 



Using the tree result of the gap equation (|37|) , we obtain 
1 



+ 



IIi 

n 2 



d D i 



(2tt) d l 2 + (m + e) 2 
i ^l-X D g ( ff(x)+4X D q D (x, S )^ 



dt 



1/2 

-1/2 
X + S 

s 





J (m + e) 2 +p 2 /A \ 


d D l 


I //4 J 



(2vr)^ [/a + p 2( 1/4 _ 4 2) + ( m + e )2] 2 

(53) 
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where 



qi(x,s) 



, 1 + x 1 + 2x + 2s , Vl + x + s + v/s 

In 1 . = In = — 

x 2ys(l + x + s) VI + x + s - Vs 

X + S . y/x~+~S + Js 

ln ; f 

s V 1 + s — v s 



(54) 



930, S) 



Va; + s(l-t 2 



: tan 



\/l + X + S + -v/s 

2y/s{l + x + s) 111 VvTTxTs- 



In 



with 



442 



< s < 



(55) 



(56) 



and Ad being given by ea.l|42|l. 

It should be noticed that in eq. the I5 part of 1/A + II2 , the (inverse) propagator 
of pions, vanishes when s 1—* as well as e 1—* 0; that is, pions are the massless Nambu- 
Goldstone particle(s). The quantity e, therefore, plays a role of an infrared cutoff. 

The onedoop part of the gap equation is derived from a restricted (onedoop) effective 
potential obtained by putting U n = in ea.(|44|). 



V ilzir=o 
where 



A D 



/•1/4 

/ dss^^ 2 ^ 2 Qi)(x, s) + ^independent terms , (57) 
J 



2iV 7r^/ 2 r(L>/2) 

Qi(x,s) = ln ( 1 - \4gf\x) + 4A 4 (x + s)g 4 (x, s 



+ ln(l - A 4 54 0) (x) + 4A 4 s<? 4 (x,s)) , 

Q 3 0,s) = In (l - A 3 ^ 0) (x) +4A 3 (x + s)q 3 (x,s) 
+2ln(l-\ 3 gl >) (x)+4\ 3 sq 3 (x, S ) 



(58) 



(59) 



to give 



Ox 



A D 



XL-=0,J=0 



2Ntt d / 2 T(D/2) j 
2N 4 d / 2 - 1 tt d / 2 T(D/2) 9d {X 



1/4 dss (D-2)/2?QD(x,s) 



dx 



(60) 
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A4#4°] + 4A 4 (g 4 + (x + s)q4 >x ) ^ -A 4 c^ + 4A 4 sg 4j:E \ 

1 - \ A gf ] + 4A 4 (a; + s)g 4 1 - X 4 gf + 4A 4 sg 4 / 

-5 4 °] + %4 + (x + s)q 4 , x ) -gfl + ^sq 4 ,x \ 

+ -7m — : . ( 61 ) 



eg^ + 4( v / x -e)(x + s)g 4 + 4( v / i - e)sg 4/ 

(1),, = _ 2 f 1/4 dsV - s ( A:^:; 1 ; • lA :{ (, /:{ • S.r • *),/,.,,) + 2 -A 3 4 l+4A 3S g3,/ 
-A) S \ 1-A 35 f +4A 3 (x + % 3 1-A 35 f +4A 3 sg 3/ 

r , r ^ -gg +4(g3 + (x + g )g3,,) + 4^,, \ 

-2(Vx - e) / dsVs —-- + 2—- — '■ , (62) 

A) V £ 5 3 +4(Vi -e)(x + s)q 3 eg£> + 4(^/i - e)sq 3 J 

with 



- d 1 D 

qn,x = • (63) 

In eas. (|61|) and (jHU), we have replaced A^ by the tree value, {-\fx — e)/(^/xg^ (x)) in the 

first expressions. In Figures ^ an d El we show the shape of (x) as well as g^ (x) for 
< x = (m + e) 2 / A 2 < 1. We choose two cases for the infrared cutoff e = e/A as 0.05 
and 0.1. 




Or 0.2 0.4 0.6 0.8 1 e 2 0.2 0.4 0.6 0.8 1 



Figure 1: g\ (a;) (a solid line) and g\ (x)(a dotted line) are plotted with e being set 0.05(left) 
and O.l(right). x = (m + e) 2 / A 2 and e = e/A. Note that gf\x) is positive but g± (x) is 
negative everywhere. 

The gap equation, up to the one-loop order of the auxiliary fields, therefore, is found 

whose right hand side diverges at ^fx = e, implying that Xd = 0. As was discussed before, 
this is physically reasonable since we have a current quark mass e even in a free A = 
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Figure 2: g% (x)(a solid line) and g\ (x)(a dotted line) are plotted with e being set 0.05 (left) 
and O.l(right). Note that g$\x) has a zero at x = 0.241(e = 0.05) or at x = 0.276(e = 0.1). 

theory. In view of ea.(|64|). an important role of the current quark mass e under the loop 
expansion is recognized: when e = (e = 0) the right hand side reads 

( 5 g ) (-)+^ i) (-))L_ o , 

whose second term becomes infinite. (Note that in eqs. (|61|) and (|62j). qo,x is singular at 
x = under e = 0.) Thus the second term surpasses the first, causing a breakdown of 
the loop expansion. On the contrary, if e 7^ 0, the second term is much smaller in the 
dangerous region ^fx ~ e owing to the factor yfx — e in front of the integrals in eqs. ljfiljl 
and (fH2f). 
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Figure 3: The ratio (g^\x)/N)/g^\x): the left graphs are for e = 0.05 and the right for 0.1. 
The dotted line stands for N — 1 and the dash-dotted line for N = 5. It is recognized that the 
ratio remains less than unity even in N — 1. 
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Figure 4: The ratio (g^\x) /N) / g^\x): the left graphs are for e = 0.05 and the right for 0.1. 
The dotted line stands for N — 1 and the dash-dotted line for iV = 5. It is recognized that the 
ratio remains less than unity even in N — 1. 



In Figures El and ^ we plot the ratios (g D /N)/g D ' in D = 4 and D = 3, respectively. 
It is recognized that even in N = l(the dotted line) the loop expansion is legitimated; 
since the ratio remains less than unity. It is also shown that the smaller e becomes the 
greater the ratio goes. The critical values that cause the ratio to exceed unity under 
N = 1 are then found such that e = 0.0326 at x = 1 in D = 4 and e = 0.040 at x = 0.0052 
in D = 3. Those of e can be put smaller when N goes larger. 

Finally we plot the right hand side of the gap equation (J64|) for D = 4 and D = 3 in the 
Figures 03 and H3 The horizontal line again stands for the value of x = (m + e) 2 / A 2 and the 
vertical line for 1/Xd- The solid line designates that of y / xg^(x)/(- v /x — e), namely the 
tree or N = oo case. The dotted and the dash-dotted lines include 0{1/N) contribution 
with N = 1 and N = 5, respectively. We have set e = 0.05(left graphs) and e = 0.1 (right 
graphs). In D = 4 magnified figures for < x < 0.15 are shown together with the whole 
plots inside those. 

It is seen that for a fixed four-Fermi coupling Ad, that is, with respect to a (supposed) 
horizontal line, mass x is a monotone increasing function of ./V in D = 4, but in D = 3 
the dependence is not so simple because of zeros in g<p (x) (see the Figures ^ and [2| ; x is 
monotone decreasing(increasing) in a small(large) mass or a four-Fermi coupling region. 
Physically speaking, due to quantum effects, xSB is restored in D = 4 at any coupling. 
Meanwhile in D = 3 it is restored (enhanced) in a strong(weak) coupling region. 



4 Discussion 

In this paper we have examined the higher order(= quantum) effect of auxiliary fields 
to the gap equation in the NJL model. Contrary to the observation by Kleinert and 
Bosschejllj. we find that auxiliary fields still play a significant role under the nonvanishing 
current quark mass e. "Pions can still survive" in the NJL model. We cannot put the 
intrinsic fermion mass into zero but at the order of 4/100 when N = 1, in order to 
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Figure 5: The right hand side of the gap equation in D = 4: the solid line designates the tree 
order or N = oo , while the dotted and the dash-dotted ones include the one-loop effects with 
N — 1 and N = 5, respectively, e is put 0.05(left) and O.l(right). Graphs with < x < 0.15 
are shown inside which the whole shapes < x < 1 are given, x is recognized as a monotone 
increasing function of iV with respect to a (supposed) horizontal line, namely a fixed four-Fermi 
coupling. 



overcome the infrared divergences and to ensure the loop expansion. In D = 4 the 
dynamical mass re is a monotone increasing function of N at any fixed four-Fermi coupling 
constant. But in D = 3 x is monotone decreasing(increasing) under a small(large) mass or 
coupling region. In other words, dynamical mass shrinks by means of quantum effects in 
the strong coupling regime of D = 3 as well as in D = 4. On the contrary it swells under 
the weak coupling regime in D = 3. We have already encountered a similar situation in 
D = 3 in the reference t 5j dynamical mass is a complicated function of the magnitude 
of the background magnetic fields(MBMF) under the influence of quantum gluons. In 
a small mass or coupling region, it is a monotone increasing function of MBMF, but 
decreasing function in a larger mass or coupling region. 

In this way, we recognize that the auxiliary field method for the NJL model can survive 
with an infrared cutoff. The power of the auxiliary field method is shown in ref. jlOj - 
using 0- and 1-dimensional examples. Cases have been, however, only for bose ones 
so that analysis for fermionic models is necessary. The 0-dimensional fermionic model, 
the Grassmann integration model, is studied to fulfill our expectation that inclusion of 
higher-loop effects of auxiliary fields improves a result much better [T7|. The 1-dimensional, 
quantum mechanical case is now under study. 
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Figure 6: The right hand side of the gap equation in D = 3: the solid line designates the tree 
order or N = oo , while the dotted and the dash-dotted ones include the one-loop effects with 
N — 1 and N = 5 respectively, e is put 0.05(left) and O.l(right). In x < 0.241(0.276) a; is a 
monotone decreasing function of N for a fixed coupling, but on the contrary in x > 0.241(0.276) 
that is an increasing function. 
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